AEI-2004-047 



One loop photon-graviton mixing 
in an electromagnetic field: Part 1 



Fiorenzo Bastianelli 

Dipartimento di Fisica, Universitd di Bologna and INFN, Sezione di 

Bologna 

Via Irnerio 46, 1-40126 Bologna, Italy 



Christian Schubert 

Max- Planck- Institut fiir Gravitationsphysik, Albert- Einstein- Institut, 
Muhlenberg 1, D- 144^6 Potsdam, Germany 

and 

Department of Physics and Geology 
University of Texas Pan American 
Edinburg, TX 78541-2999, USA 



Abstract 

Photon-graviton mixing in an electromagnetic field is a pro- 
cess of potential interest for cosmology and astrophysics. At the 
tree level it has been studied by many authors. We consider 
the one-loop contribution to this amplitude involving a charged 
spin or spin 1 /2 particle in the loop and an arbitrary constant 
field. In the first part of this article, the worldline formalism is 
used to obtain a compact two-parameter integral representation 
for this amplitude, valid for arbitrary photon energies and back- 
ground field strengths. The calculation is manifestly covariant 
througout. 



1 Introduction 



In the presence of an external electromagnetic field, many quantum pro- 
cesses exist which are forbidden in vacuum (see, e.g., 012]). In particular, 
transitions between bosons of different spin become possible [3]. One 
such process which has been well-studied is the axion-photon mixing in a 
magnetic field [SIEIIZI- Similarly, also photon-graviton mixing is possible 
in an external field [HH^E]. The corresponding tree level amplitude is con- 
tained in the coupling h^^T^^ of the graviton h^y to the energy-momentum 
tensor T^^ of the electromagnetic field, 

rp^^<^ ^ F'^^F'^ - 1f„^F"V- (1-1) 

Namely, when taking F^^ = F^^^ + f^'^ with the external field strength 
tensor and /^'^ the photon field, hnuT'^'^ yields the trilinear term 




Due to the smallness of the gravitational coupling constant k this process 
has attracted less attention than the axion-photon coupling. Nevertheless, 
its relevance for astrophysics has been scrutinized by a number of authors 
(see [H] for a discussion of possible laboratory experiments). In ^ photon- 
graviton conversion near a pulsar was studied but the transition rate was 
found to be very small. In ^] it has been suggested that the same 
conversion due to a primordial magnetic field could be responsible for the 
observed anisotropy of the cosmic microwave background. However, find 
that its effect becomes negligible for standard cosmological magnetic fields if 
plasma effects are taken into account. Renewed interest in this amplitude has 
been generated by the recent models with large extra dimensions 12 . These 
models contain additional massive Kaluza-Klein gravitons which might lead 
to an enhancement of the photon-graviton conversion effect. See |13j for 
a discussion of possible observable effects in astrophysics as well as in the 
laboratory. In the same context also graviton-photon conversion on spin 
and 1/2 particles was considered jl4j . 

To our knowledge, the photon-graviton amplitude has so far been consid- 
ered only at tree level. In the present paper, we extend its study to include 
the one- loop corrections due to virtual spin or 1/2 particles. Contrary to 
the tree level case, at one loop these amplitudes depend nontrivially on both 
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the photon energy and the background field strength. Thus it is conceivable 
that, for some region in this two-parameter space, the one-loop amplitudes 
might be comparable with or dominating over the tree level ones. In the first 
part of this article, we use the 'string-inspired' worldline formalism to ob- 
tain compact integral representations for these amplitudes. Their numerical 
study will be undertaken in the second part. 

Our calculation provides also a first example for the application of the 
string-inspired technique to mixed photon-graviton amplitudes. In recent 
years, methods derived from or inspired by JS] string theory have been 
extensively used for the calculation of on-shell gluon and graviton ^1 
[T9] amplitudes, as well as for QED amplitudes in a constant field [201 
m m m HSl Uni I2Z1- More recently, in |^ also the contribution to the 
graviton vacuum polarization involving a scalar loop was obtained. The 
treatment of gravitational backgrounds in the wordline formalism involves 
a number of mathematical subtleties which have been clarified only recently 

US EOl ED ElESl El • 

Phenomenologically, the photon-graviton amplitude is primarily of inter- 
est in the magnetic field case. However, in the formalism used here it makes 
technically no essential difference whether one calculates an amplitude in 
a constant magnetic or in a general constant field. Thus we will keep the 
electric component. 

The organization of the paper is simple: In chapter |2l we present the 
general formalism for calculating one loop amplitudes involving either a 
scalar or spinor loop and any number of photons and gravitons, in vacuum or 
in a constant external field. Chapter El contains the scalar loop calculation, 
chapter H the spinor loop one. Our conclusions are given in El In the 
appendix we verify that the results obey the relevant gauge and gravitational 
Ward identities. 



2 Mixed electromagnetic - gravitational amplitudes 
in the worldline formalism 

The application of the worldline formalism to flat space calculations has 
been described in detail in the review jSH]. However, the generalization 
to processes involving curved backgrounds, or gravitons, is less familiar. 
Therefore, after a brief description of the worldline formalism adapted to 
the case of a constant electromagnetic background, we will describe certain 
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subtleties arising in this formalism with the coupling to gravity. 

Let us first consider the case of a scalar particle coupled to electromag- 
netism and gravity. We use here euclidean conventions. The QFT for this 
scalar particle is described by a complex scalar field (j) with action 



(2.1) 



where ^ parameterizes an additional non-minimal coupling to the scalar 
curvature^. The corresponding one- loop effective action is formally given by 
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T[g,A] = - \ogT)ei-^ {-U A + iR + m^) = Ty:\og{-n A + iR + m^) 

(2.2) 

where VIa is the gauge and gravitationally covariant Laplacian for scalar 
fields obtained from 1)2. This effective action can be represented in the 
worldline formalism by 

n9,A] = -r^f Pxe-^[-^^'^l (2.3) 

Jo -t JPBC 

with the worldline action 

S[x;g,A]=J^ dT[-gf,^{x)i''x'' + ieA^{x)±'' + ^R{x) + m'^j. (2.4) 

This worldline representation contains a standard integral over the proper 
time T, and a quantum mechanical path integral over the particle coordi- 
nates x^{t). Due to the trace in (|2.2)) the path integral is to be taken with 
periodic boundary conditions x'^(O) = x^{T). Thus it corresponds to an 
integral over closed loops in spacetime. In the worldline or 'string-inspired' 
formalism, the path integral ()2.3() is manipulated into Gaussian form and 
then evaluated using appropriate worldline correlators. In an amplitude 
calculation, the external fields are specialized to plane waves. Then each 

^The value ^ = is the minimal coupling, while the value ^ = ^(^^i-^ gives a confor- 
mally invariant coupling in the massless case. 

■^Our present definition of the effective action differs by a sign from the conventions in 
however, there is no difference at the amplitude level. 
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external leg is represented by a vertex operator. For the electromagnetic 
field this is the photon vertex operator 



VAAk,e] = ej^ dTx^{Ty^<^) (2.5) 

with an associated coupling constant —ie. The graviton vertex operator will 
be derived below. For the calculation of the path integral, one first splits off 
the loop average position: 



x'^(t) = x^ + y^^), (2.6) 
^ ^£drx^{r). (2.7) 

The path integral then factors into J Dx{t) = J d^xo J Dy{T). The inte- 
gral over xq just produces the global delta function for energy- momentum 
conservation. The reduced path integral / Dy^r) is Gaussian and can be 
calculated using the Wick contraction rule JS] 



{y^iri)y''iT2)) = -5'^"GB(Ti,r2). (2.^ 
Here is the 'worldline Green's function' 



^ / N I I (n - r2) 



T 

Ctb{ti,T2) = sign(ri - T2) - 2 

GB{ri,T2) = 26{n-T2)-^. 



(2.9) 



Here and in the following we will often abbreviate Gbi2 = Gb{ti-,T2) etc., 
and a 'dot' generally refers to a derivative in the first variable. 

The inclusion of a constant electromagnetic background field F^^ can 
be achieved easily [201 using Fock-Schwinger gauge centered at the loop 
average position xq. In this gauge 



A^,{x) = ^y^F,^ (2.10) 
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so that the presence of the background field produces only an additional 
term 



AL = hey'^F^^r (2.11) 

to the worldline Lagrangian in (|2.4() . Since this terms involves y only 
quadratically it can be taken into account by an appropriate change of the 
worldline correlators. Namely, instead of (|2.(S|2.9j) one finds 

{y^iTi)y''iT2)) = -ar(n,r2) (2.12) 
with a field dependent worldline correlator |37| I23j 



Z \sin(2^) 

gB{ri,T2) = 25i2-|-^e-*^^^- 

1 sm(Z) 



(2.13) 



where Zf^u = eTF^y. These expressions should be understood as power 
series in the Lorentz matrix Z\ more explicit expressions are given below. 
Note the symmetry properties 



Qbi2 — 0b2i^ Qbi2 — —Gmi^ Qbi2 — Gmi- 



(2.14) 



Contrary to the vacuum case, in a constant field background the worldline 
correlators have non-vanishing coincidence limits. 



Ob{t,t) = -L_(^z cot{Z)-l 
Gb{t,t) = icot{Z)- — . 



(2.15) 
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In the following a 'bar' on a quantity denotes the subtraction of its coinci- 
dence limit, e.g., 



Gbi2 = ^512 - Gbu ■ (2.16) 

Note that coincidence limits of worldline correlators are always constant due 
to their translational invariance. 

The path integral determinant also becomes field dependent. In flat 
space ,20^ 



j -Dy(r)exp - dr^^ 



sin(2^) 



Z 
(2.17) 

The path integral representing a spin ^ particle in an electromagnetic 
field differs from the spin case 1)2. 3|) above by a global factor of — ^, and 
by an additional Grassmann path integral representing the spin. 



ABC 



Vil) exp 



(2.18) 



Here the path integral is over antiperiodic Grassmann functions, tp^^{T) 
—'ijj^{0). In the vacuum case, the appropriate worldline correlator is 



(V'(ri)V'(r2)) = ^GFiTi,T2) = ^sign(ri - T2) . 
For a constant background field this correlator turns into 



(2.19) 



{riri)rir2)) 



1 



^F''(n,'r2) 



where 



(2.20) 



GFin,r2) 
Qf{ti,T2) 



-iZGi 



Gp 



12- 



cos(Z) 



Gfi2 + Gp 



2iZ 



12 



Tcos{Z) 



-iZGBi2 



(2.21) 
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Its symmetry properties are 



Gf{ti,T2) = -^F(T2,ri), Gf{ti,T2) = Q'§'{t2,Ti) . 



(2.22) 



The coincidence limits are 



Qfii = — itan(Z) , Qfu = 2(^ii + 2eFtan(Z) 



(2.23) 



The free path integral in a constant F^y background in D dimensions is 
normalized as 



j Dipexp — J dT(^—'ip-ip 



o , 1 

2Tdet2 



cos(^) 



(2.24) 



The photon vertex operator (|2.5|) acquires an additional Grassmann piece, 



dr 



x"(r) + 2i'0"(r)'0(T) • k 



^ik-x(T) 



(2.25) 



In the flat space case, the naive Gaussian path integration gives unam- 
biguous and well-defined parameter integral representations. To the con- 
trary, when the coupling to gravity is introduced more precise definitions of 
the worldline regularizations are required to define the path integral prop- 
erly. Let us describe these issues briefly. First of all it is convenient to 
exponentiate the nontrivial path integral measure in a regularization inde- 
pendent way by using ghost fields [2^1 OOI • The covariant measure in 1)2. 3|) 
is of the form 



Vx = Dx ^deigf,i,{x{T)) 



(2.26) 



0<r<T 



where Dx = Hr dPx{T) is the standard translationally invariant measure. 
It can be represented more conveniently by introducing commuting a'^ and 
anticommuting b^,c^ ghosts with periodic boundary conditions 

Vx = Dx TT Jdetg.^{x{T)) = Dx f Z)aL>6L»c e-^s'^!^''''*'^! (2.27) 

0<r<l ^ -^PBC 
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where the ghost action is given by 
Sgh[x,a,b,c\ 







(2.28) 



The extra vertices arising from the ghost action guarantee that the final 
result will be finite (see, e.g., f5Bj). 

For the perturbative calculation of graviton amplitudes around flat space, 
one next linearizes the metric. 



and then specializes h^^{x) to plane wave form. 



(2.29) 



h^v{x) 



ik-x 



(2.30) 



This leads to the following vertex operator for the graviton coupled to the 
loop scalar 



'flU 



dr 



x^(t)x'^(t) + a^'{T)a''{T) + 6^(T)c'^(r) 



J,k-x{T) 



(2.31) 



with an associated coupling constant factor of — -j. The Wick contraction 
rules for the ghosts are 



(a^(n)a^(T2)) 
(6^(n)c^(T2)) 



25{n - T2W\ 



(2.32) 



Similarly, for the fermion loop case one finds a graviton vertex operator 



' spin 

+2 



dT 



x^(t)x^(t) + a^'{T)a!'{T) + 6^(T)c'^(r) 



(jJ^'{T)iJ''{T) + a^'{T)a''{T) + ix^'{T)^'' {t)^{t) ■ k 



(2.33) 
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where are the additional bosonic ghosts arising from the nontrivial path 
integral measure for ip^ [23] . Their correlator is 



(a'^(ri)a-(T2)) 



(2.34) 



In the perturbative expansion around flat space various worldline Feynman 
diagrams are linearly and logarithmically UV divergent. In fact the ad- 
ditional derivative couplings due to the metric worsen power counting in 
momentum space. The contributions from the measure, i.e. terms involving 
the ghost correlators, will always eliminate these divergences. Nevertheless, 
finite ambiguities are left over and dealt with by specifying a regulariza- 
tion scheme together with renormalization conditions ^ . These amount to 
the requirement that the path integral in (|2.3() correspond precisely to the 
Hamiltonian operator H = {—Ca + + C-R) appearing in (|2.2|) . These 
renormalization conditions produce a finite counterterm of the form 



that must be added to the action (|2.4j) . Three regularization schemes have 
been worked out in detail for this purpose: mode regularization (MR), time 
slicing (TS), and dimensional regularization (DR). The corresponding coun- 
terterms are given by 



For details see jSl EH Ell- 
in the present article, we will be interested in the graviton-photon ampli- 
tude in a constant electromagnetic background mediated by a virtual scalar 
particle loop. This means that we will only need to consider the linear 
coupling to the metric fiuctuations around flat space. To this order the 

^We emphasize that these issues concern the one-dimensional worldline theory and 
bear no direct relation to the issue of regularization in spacetime, to be considered later. 




(2.35) 



Vmr 



Vdr 



Vts 




(2.36) 
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differences between the counterterms in ()2.c{6|) can be neglected since they 
are at least quadratic in metric fluctuations. On the other hand the leading 
part of the counterterm contained in the —^R piece effectively changes the 
coupling ^— >C = ?~iiii (|2.4() . This implies that to this order all three 
regularization schemes can be used interchangeably. 

In the case of the spin ^ particle the coupling to R is fixed by the Dirac 
equation and corresponds to ^ = 0. Additional details on the worldline 
formalism with background gravity can be found in [28L I35| . 



3 Calculation of the photon — graviton amplitude 
in a constant electromagnetic field: scalar loop 

According to the above, the amplitude for the interaction of a graviton h and 
a photon A via a scalar loop, in the presence of a constant electromagnetic 
background field F^^^, is given by the following expression: 



{h{ki)A{k2)) 



Here represent the photon and graviton vertex operators for the scalar 
loop case, (^3]) and (|OT|) . 

To start with, we perform the split (|2.6() and the trivial xq integration 
which produces the global 5 - function for momentum conservation. The 
remaining path integral / Dy(T) is Gaussian and thus can be reduced to 
Wick contraction. Taking (|2.17|) into account, one obtains 



^ 4' Jo T 



-m?T 



X l^DxDaDhDcVlU[ki,e'^]V,ti[k2,e' 



xexp 



/ dri-ix^ + + h- c) + -iex^'Faux'- 
jQ ^4 2 



(3.1) 



{h{ki)A{k2)) = —{2TTf5{ki + k2) — e-'"^(4^r) 2det-2 ^ ^ 
4 Jo 1 

x{v,t,[kue'']V±,[k2,e' 



(3.2) 



eTF, 



flU) 
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Performing the Wick contractions according to the rules ()2.8|) . H2.32() . 
and using ()2.14|1 . we obtain 



{h{kM{k2)) = {27:)^5{k, + k,)el,etli^^^J{k) 



(3.3) 



where k = ki, and 



(k) 



— g-m ^T--det"2 



4(47r)Tio - 

X / dTi I 
JO Jo 



sin(Z) 



A-r^p-k■QBl2■kT^f^,a 



(3.4) 



= -{G'mi--^5ii5'''){k-QBi2y-[Q%\,{QBi2-k 

+ (Gbi2 ■ ky[gBi2 -kjik- GbuY - 4C(5^"fc2 - pr)(fc • Gbu) 

(3. 

It is useful to add to Z^^'" g-fc-Gsi2-fc ^j^g total derivative term 



2 9ri 

Then /^''^X gets replaced by J^^^'f, 



(3.6) 



"^scal ~ "^scal,! "i" '^scal,2 "i" "'scal,3 "^sca 



a 

seal, 4 



where J^^'^) = Ujt^'^ + jn^) and 



(3.7) 



•^scaLl 



TfJ.u,a 
'^scal,2 

•^scal.S 



- ^b"2 ( Gb12 ■ k 

■ kY\{QBi2 -k) (k- Gbu 



—fia / .. 
QBl2\pBl2 ■ k 



'B12 



a —VOL 



'5b\2^ • ^512 • ^ 



-^i[b^-k^ -k^k''){k- Qbx^ . 



(3. 
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Before proceeding further, let us use this integral representation to analyse 
the general structure of this amplitude. Although our calculation is nonper- 
turbative in the external field, we can, of course, use the series expansions 
of the worldline Green's functions (|2.13|) to compute the amplitude involv- 
ing a given number of interactions with the field. It is then immediately 
seen that this amplitude is nonzero only if this number of interactions is 
odd, since otherwise the ti^2 integrations vanish by antisymmetry. Also, 
the integrand contains terms which are singular at T = 0, indicating UV 
divergences. While so far our calculation has been valid for any spacetime 
dimension D, the structure of these divergences depends on D. We therefore 
confine ourselves to the four dimensional case in the following. For D = 4 
the terms ^^^i" 4 finite, while Jg'^ail 2 contain terms with a loga- 

rithmic divergence at T = 0. Those divergent terms involve the field only 
linearly, and are thus easy to compute by expanding the formulas ()2.13() to 
the linear order in F. In dimensional regularisation, the result is 

where 



(3.10) 

As expected, this counterterm is just the momentum space version {f^^, = 
kfM^u — kySp) of the tree level interaction term (|1.2|) . 

We perform renormalization by subtracting the amplitude at zero field 
strength and zero momentum. This can be done under the T - integral, 
leading to the following form of the renormalized amplitude II: 





■sin(2^) ■ 




|det"^ 


Jo Jo 




^ z \ 



(3.11) 



Further, again following [^1^2] we split QbiJ into 



Gb = Sb + Ab 
12 



(3.12) 



where S'^ contains the even powers of F^^'^ in the power series representation 
of , and the odd ones. After this replacement all terms in the 
integrand are either symmetric or antisymmetric under the exchange ri <-> 
r2, and the antisymmetric ones can be deleted since their ti^2 - integrals 
vanish. Note that the exponent k ■ Q bi2 ■ k = k ■ Sbi2 " ^ is symmetric under 
this exchange. 

Further, as usual |3ni we rescale Tj = Tui, i = 1,2, and use the translation 
invariance in r to set U2 = 0. After analytic continuation to Minkowski 
spacetime, this leads to our final result for this amplitude: 



^kclW - g4^2 



det 2 



s 



■sin(Z)- 


L 


I z \ 





du-i 



m=l 



(3.13) 



where C^'^'" is as in (pnU)) (with 6'"' rji"') and 



'^scal,2 



seal, 3 



'-'B12 



(Abi 



—na 

+ A 



12 • k 

512 ('^512 • k 



AWl2{SB12-k 



Sbi2 • k 



- S'^'12[Abi2 ■ k 
Sbi2 ■ k) (k- Abi2 



+ 



[Ab12 



k] (k- Sb 



12 



-•^312^ ■ Sb12 ■ k 



+ 



[Ab12 



■ k 



•^scal^ 



5^12 ■ kj l^k ■ Sbi2 
-S^Bl2k-SBi2-k 
A^{'n^"'k^ -k^k'')(k- Jbi2 



+ [ABi2-k) (k- Ab 



312 



(3.14) 



For applications of this amplitude it will be necessary to write the matrix 
functions S, A and their derivatives in more explicit form. A suitable rep- 
resentation has been given in [IHlEni- Let = ^(S^ - E"^), ^ = B • E the 
two Maxwell invariants and 



We use the metric tensor (j?^") — diag( — h++) in Minkowski spacetime. The analytic 
continuation from Euclidean to Minkowski space amounts to substituting S'^" — > r;*"^, 
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so that - 6^ = - E'^ and a6 = E • B. Let 



(3.15) 



Z-y. = iesa, = —esh. 

Then the determinant factor decomposes as 



(3.16) 



Let further 



det 2 



sin(2^) 



sinh(z+)sinh(2:_) 



aF-bF 
a2 + 62 ' 



Z_ = -i 



bF + aF 
o^ + b"^ ' 



(3.17) 



(3.18) 



with F^y = 5£^i/a/3-P"^ the dual field strength tensor^. The matrices Z± 



fulfill 



and 



Z+-Z. 



(3.19) 



F2 - 



Z2 =- 



F^ + aH 
a? + b'^ 



a^ + b^ ' 

Then one has the following orthogonal matrix decompositions: 



(3.20) 



Sb12 
•5312 



a=± " 



a=± 

"4b12 = X] ^512-2^0' 

a=± 



^In our conventions 
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<Sb12 = Gbu^ ^a^B12^a) 

* a=± 

-2 

Ab12 = - V ZaS%i2Za. 

a=± 



(3.21) 



These formulas are written in terms of the following four basic scalar, di- 
mensionless coefficient functions: 



'-'B12 



sinh(z±GBi2) 



sinh(2;±) ' 

_ cosh(z±G'Bi2) _ J_ 
~ sinh(z±) z± 

(3.22) 

Note that with respect to the exchange u\ U2 the S^^2 ^^^^ O'^'^ ^'^'^ 
^Bi2 6ven. Thus only the latter have non-vanishing coincidence limits, 

= coth(z±)- — . (3.23) 



4 Calculation of the photon — graviton amplitude 
in a constant electromagnetic field: spinor loop 

The corresponding calculation for the spinor loop case proceeds in a com- 
pletely analogous way: 

{h{k{)A{k2)) = --{-ie){--) r—e-"^^^ [ DxDaDbDc [ Dil^Da 
2 4 Jo T Jp J A 

xF^jA:i,£'^]yAjA:2,e^ 

r rT 

xexp 



- / c?r(^(i;2 + + 6 • c) + hex^'Fr.^x" + ^(V' • V' + a^) - ie^l^'^Fau^ 
JO ^4 z Z 

= — — 2 2-(27r)^(5(A;i + A;2) / — e'"" ^ (AttT) 2det"2 ^ ' 
8 Jo T 



Z 



(4.1) 



15 



Here V^l^ now represent the photon and graviton vertex operators for the 
fermfon loop case, (|2.25|) and (|2.33|) . The additional Wick contraction rules 
have been given in H2.20() . (|2.34() . Performing the Wick contractions in H4.1|) 
we obtain the analogue of H3.4() . 



Here 



where 



D_ 

eK2 2 



-^^T-fdet-H*^ 



8(47r)--^o I 
Jo Jo 



Z 



spin 



spin 



0) + jC/^'^)." 

^) ^ -'extra 



(4.2) 



(4.3) 



extra 



+ 



[QF12-kYgTl2-GTl2{QFl2-k 



+ 



+ 



+(g^\^-25i^5^'^)[k-gBi2y 

2hi6^''' -(GBi2-kY(gBi2-kY][ Qf22 ■ k 

Qfw ■ k 



^Bll 



Gbi2 ■ k]^( k - Qb12 



(4.4) 



It is useful to replace by 



Cr-^srar(e = o) + /, 



extra 



(4.5) 



The integrand can then be rearranged in the following way (compare 
with ()3.8() for the scalar loop): 



TfJ.u,a 
spin 



spin,l ' spin, 2 ' spin, 3 



(4.6) 
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where 



spin,l 



•^spin,2 



spin, 3 



{^Bll ~ ^Fll) {Qb21 + Qf22) ■ k 

Gbi2(Sbi2 ■ - ^^12(^^12 • fc) - Gbi2 {0bi2 + Gfu) ' k 

+^F12 {Qf12 ■ k^ , 

-{GBi2-kj ^^{Gbi2 + Gfu) ■ k^ |^(^B21 + ^^22) • A; 

QFi2kj \QF2ikj +GBi2k ■ Gbi2 ■ k - gp°l2k ■ Gfi2 ■ kj. 

(4.7) 



As in the scalar loop case, renormalization requires only the subtraction of 
a logarithmic divergence, yielding 



n^pin (k) 



en 



dT 



X <^ det 2 



327r2 Jq r3 
tan(^)^ 



Z 



T A 



(4.8) 



The tensor C^^'" appearing in the counterterm is the same as in the scalar 
loop case, eq. 

We now need also the analogue of the decomposition formulas 1)3. 12(1 . 
(Uni), 'l3:22]> for Gf, 



Gf 



(4.9) 



5^12 
•Af12 

Sf12 

•Af12 



" X] ^F12 1 
a=± 

-i ^ ^F12 

a=± 

2i 



Gf12 1 Y] ZaA%i2 ^1 , 

2 



a=± 



F\2 J 



(4.10) 



a=± 
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written in terms of the basic coefficient functions 



'-'F12 



= Gf 



12 



= Gf\2 



cos\i{z±Gbi2) 
cosh(z-i-) 

sinh(z±G'Bi2) 



cosh(2;±) 
tanh(z±) . 

As in the bosonic case S'^^^^ (^fi2) ^'^'^ (even) with respect to 1 2. 

Proceeding as in the scalar loop case, we obtain our final result for the 
amplitude: 



eK 



ds 



jfj,v,a 
spin,3 



327r2 Jo s 



tanh(z+) tanh(2;_) Jo 



du^ 



m=l 



(4.11) 



where 



spin,l 
'^spin,2 



- ('^^ll - <5^1l) {k ■ (^512 + ^F22))" , 



-AWi2 ( Sbi2 ■ky + AWi2 (5fi2 • k) ' 

-^312 (<5b12 • k^ - A^^2 {Sf12 ' k^ 

S^Bu ( Ubi2 + ^Fii) ■ky + S'f'I^ (-4fi2 ■ kf , 
SBi2-kY (^SBi2-k^ (^k- {Abi2 + Afii)) -(^SFi2-k^ ( 
+ {{Abi2+Afii) -kyik- SbuY - (Afi2 -kjik- 5^12)" 

—va 

- ABuk ■ Sbi2 ■ k + AF{2k ■ Sfi2 

+ [Ab 



+ 



[k ■ Afi2) 



/iB12«; ■ OB12 ■ K -r ^F12'^ ■ ^F12 " K 

^ ^ Lv'^-^^^ " ^) ' '^^12) - (Sfi2 -k^ [k- SF12) 
■{{Abi2 + Afu) ■ kY (k ■ {Abi2 + AFii)y - (AFi2-kY{k- ^^12)" 
-S^'l2k ■ i'-L -c'^" u. c u] 



■ 5si2 • k + Sp^2k ■ Sfi2 ■ k 



(4.12) 
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5 Conclusions 



We have obtained compact explicit integral representations for the one-loop 
photon-graviton amplitudes involving a scalar or spinor loop and a constant 
electromagnetic field. The use of the constant field worldline formalism along 
the lines of [SSj has allowed us to achieve this with modest calculational 
effort, and without having to specialize to a special Lorentz frame. As 
usual in this formalism the calculation for the spinor loop case has been an 
extension of the scalar loop one. Since the application of this formalism to 
gravitational backgrounds is still novel we have presented the calculation in 
some detail. 

We have verified that our results for these amplitudes obey the gravita- 
tional and gauge Ward identities. As a further check, we have also shown 
that their low energy limits agree with the result of an effective action cal- 
culation 

In the upcoming second part of this paper we use our results for a nu- 
merical study of the photon-graviton conversion process in a magnetic field. 

Acknowledgements: C.S. thanks H. Gies for helpful discussions. Both 
authors thank Stefan Theisen and the Albert-Einstein Institute, Potsdam, 
for hospitality. C.S. also thanks V. Villanueva and the Institute of Physics 
and Mathematics of UMSNH for hospitality. 

6 Appendix: Ward identities 

In this appendix we derive the relevant Ward identities and verify that our 
graviton-photon amplitudes satisfy them. This provides a good check on 
the correctness of our calculations. 

There are two types of Ward identities: one that originates from gauge 
invariance and one that follows from reparametrization invariance (general 
coordinate invariance) . 

Gauge transformations are defined by 



with an arbitrary local parameter A. Then gauge invariance of the effective 
action 




6Gr[g,A]=0 
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implies that 



/ 1 6T 



. 



(A.3) 



Similarly, infinitesimal reparametrizations are given by 

6rA^ = ed^A^ + d^eA, , 5Rg^, = V^Cu + V.^/. (A.4) 
with arbitrary local parameters The invariance of the effective action 

6RT[g,A]=0 (A.5) 

now implies 

1 6T 



^\^/9 5g^lu ^/g 5A^ 



-VA,, 



. 



(A.6) 



The Ward identities thus obtained can be combined and written more con- 
veniently using standard tensor calculus as follows 



5T 

saT,, 



0. 



(A.7) 



2d, 



6r 



+ 



6r 
Ia:, 



d,A- + r-,i2 



6r 



— -0 



Now we are ready to consider the special case of the graviton-photon cor- 
relation function in flat space and in a constant electromagnetic background 
Ff^iy described by the gauge potential ^/^(x) = ^x'^Fi,^ 



(A. 



Taking functional derivatives on the general Ward identities ()A.7|) to relate 
them to eq. (|A.8|) we obtain 



(A 



5^T 



^ ^-'V)^ 5A^{x)5A^{y) 



5^ 

F ^ A 

^ 5A^{x) 



Now we Fourier transform these identities to momentum space 



j dxi..dxn e 



ik-LXi+..+ik„Xn ■ 



{xi,..,Xn) 



(fcl,..,fe„) 
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and from eqs. ()A.9|) we obtain 



fe"rfr^) = o, (A.io) 

One may already notice that the term proportional to '^^q-^ can be discarded, 
since it vanishes at zero momentum. Thus we see that the gravitational 
Ward identity relates the graviton-photon amplitude to the photon-photon 
amplitude. The latter has been calculated in the worldline formalism in 
PSI IHH] . The two point functions used above are simply related to the 
vacuum polarizations computed in the main text and in |25| 136) as follows 

rf^"_^^ = -n'^'"(A:), 

«rf;;^,) = -m'^'^ik). (A.ll) 

Thus the expected Ward identities are 



(A.12) 



and 



-nF\Il^%k). 



(A.13) 



Let us start with the amplitude due to a scalar loop. These Ward iden- 
tities are most easily checked using the form ()3.8p . To verify the gauge 
Ward identity ()A.12|) note that 3 4 vanish at the integrand level when 

contracted with ka. This is not the case for ^o</^ai"i' ^'^^^ term the 

integral vanishes because k ■ Qbu • k is antisymmetric in ri^2- 

The verification of the gravitational Ward identity HA.13|) requires a bit 
more work, since here only k^J^^J^^ drops out at the integrand level. We 
can simplify it by adding a suitable total derivative term: 



, jfiv,a (.) , jixv,a (.) _ / Qva _ 9 X xua\ _^ {•) 

+ -^< \k-QBVl) [GB12-kj - GbU^ ■ Gb12 ■ k 

= {(k'^Bi2y{0Bi2 + 26ul) -kj - {pBi2 + 25iityk ■ Gbu ■ k\ e^'^ 
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(A.14) 

Here in the last step we have used the fact that Qbi2 and k ■ Gbi2 ■ k have 
vanishing coincident Umits to replace 26ii by Gbu — Gbi2 = 2(5ii — 25i2- 

Next, let us write down the representation corresponding to H3.4() for the 
photon - photon amplitude in a constant field [2^1 EEI : 



(47r)T Jo 



— e^m ^r~-det"2 
T 



sin(Z) 



T 



(A.15) 



where 



-^s^cal = ^B12 k ■ Gb12 ' k - [G B12 ' k) [k ■ Q B12 



(A.16) 



Using ()A.14j) . and (|A.15|1 . the gravitational Ward identity (|A.13|) can 

now be easily verified using the matrix identity 



^512 — Gb 



12 



2ieF ■ Gb 



12 



(A.17) 



Let us now consider the amplitude due to a fermion loop. The gauge 
Ward identity is again easily seen to be satisfied. To check the more subtle 
gravitational Ward identities we need the photon-photon amplitude in the 
constant electromagnetic background |25| l36] 



where 



o / — e"*" ^r--det-2 

2(47r)--'o T 



tan(^) 



Z 











(A.18) 



spin seal extra 



(A.19) 
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with /^"[ given in ()A.16|) and 

+ {Gfi2 -ky^k- GfuY + {Gfu ■ ky{gF22 ' kY 

- {Qfii ■ky(k- Gbu) " + (Gbi2 ■ k) {gF22 ' k)" . (A.20) 

We essentially proved earlier that the scalar part satisfies the gravita- 
tional Ward identity. Thus we contract the remaining terms with k^^ and 
add a suitable total derivative 



, J{^lu),a (.) _d_rpua 
"'A'-' extra ^ ^ 



(A.21) 



with 



(^Fii • k^ {Qf22 ■ k^ + (Gfu -k^ {k- QB12) 
+ QT12 {k ■ Gfu ■ k) - {Gfu ■ ky {k ■ ^^22) " 



.(■) 



(A.22) 



to obtain 



, J(^lu),a (.) _^rj.U 

'^'A*-' extra ' ^ 



Now using 



k ■ Gf 



Gf22 ■ k 



{k-GB 

+GTi2{k ■ Gfu ■ k) - {Gfi2 ■ ky (k ■ Gf^Y 
(k-GFiiY -26uk'' (^"^ 



(A.23) 



^_Bi2 = Gbi2 + 2ieFGBi2 
Gfi2 = 2(5i2 + 2ieFGFi2 



(A.24) 



one can show that eq. HA.13|1 holds. 
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